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(54) MATHEMATICAL RELATION IDENTIFICATION APPARATUS AND METHOD 



(57) An apparatus for identifying a mathematical 
relation between base variables (x1, x2, ... f xp), when a 
set of input data d is composed of p base variables (xl , 
x2, ...,xp), and a plurality of data set d(i) of such a data 
set d is inputted, the data sets d(i) are distinguished by 
an input data distinguishing parameter i; 

The principle of the present invention is that when 
vectors in a q dimensional space, mapped from the 
input data through q base functions (f1, ...,fq), form a 
plane, a mathematical relation can is a linear combina- 
tion of the base functions; 

A set of base functions (f 1 , ... fq) are previously pre- 
pared; Then sets of the values F(i) = (F1i, ...,Fqi) of the 
base functions corresponding to the input data d(i) are 
acquired; The sets F(i) are considered as points in a q 
dimensional space; Direction cosine of a mapping plane 
is acquired through cofactors of a determinant of these 
points, or by solving an eigenvalue problem for deter- 
mining a plane, the square sum of the perpendicular 
lines from these points to the plane is minimum; When 
the direction cosine of the plane is (L1, ...,Lq), the fol- 
lowing mathematical relation is outputted: 



£(Lk-fk) = 0 
k=1 



The present invention can be used in interpolation 
of data, identification of a trajectory or form. 



Fig. I 
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Description 

(Technical Field) 

5 [0001] The present invention relates to an apparatus for identifying a mathematical relation between variables of 
measured data. An apparatus for identifying a mathematical relation identifies a mathematical relation between numer- 
ical data of, for example, a form, trajectory, etc. The extrapolation, interpolation, storing of such data can be facilitated. 
[0002] For example, the apparatus can be used as follows: 

w 1) Expression of a shape of an object: 

In a designing of a shape of an automobile, data of a shape of an object can be stored as a set of numerical data. 
However, when a shape can be expressed as a mathematical expression including a plurality of parameters, data 
can be easily extrapolated and interpolated. And new shape can be easily obtained, by varying the value of some 
parameters. And a computer analysis of the new shape, for example, aerodynamical resistance, etc, becomes 

is easy. Further, a preparation of a mold for fabrication of an object having the corresponding shape, using a numer- 
ically controlled working machine, becomes easy. 

In a shoe shop, even when a client found a pair of shoes, the size of which is identical to the size of his feet, 
and the design of which is acceptable, if the pair of shoes do not fit actually to his feet, he will not buy the pair of 
shoes. There are many variety in the shape of foot of human being. Therefore, sizes including small number of 

20 measured values of a foot are not sufficient to express a shape of a foot of human being. No apparatus or method 
is known, which can determine values of parameters of a mathematical expression, representing a complex shape 
of an object, for example, a shape of shoe or a shape of foot of human being. 

2) Identification of a mathematical relation between variables of experimental data, when the mathematical relation 
25 shall be expressed by a non-linear function: 

Even when a mathematical relation between obtained data is non-linear, data are often extrapolated and inter- 
polated, using a linear approximation. However, when a mathematical relation between the variables of data can 
be identified, data can be exactly extrapolated and interpolated, without using linear approximation. 

Assuming that a food is produced from stuffs A, B. And a characteristic value of the food X, for example, con- 
30 centration Z of an amino acid, is a function of the concentration X, Y of the stuffs A, B, the pressure P and temper- 
ature T of a treatment. When the mathematical relation between these variables can be identified, it becomes easy 
to obtain the most preferable values of the X, Y, P, T, for getting the most appropriate concentration Z, using the 
mathematical relation. Thus the development of a new product can be facilitated. 

For example, assuming that the relation between the variables X, Y, P, T, Z is as follow, no apparatus and 
35 method is known, which can determine the coefficients K0, K1 , K2, K3, by small times of experiments. 

Z/(X-Y) = K0 + K1 -P1 +K2-T + K3«X 

3) Representation of a trajectory of a moving object: For example, it is said that the trajectory of golf club head of a 
40 professional golfer is approximately a plane, on the other hand, the trajectory of golf club head of an ordinal ama- 
teur is not a plane. The trajectory of golf club head is not a line, thus the mathematical expression of the trajectory 
of golf club head is not simple. An apparatus for identifying the mathematical relation identifies the mathematical 
relation of a trajectory of such a moving object so as to analyze the trajectory. 

Trajectories of an airplane flying around over an airport, a ship sailing on the ocean, or a car running on the 
45 road, are not linear, but the changing rate of the direction of the movement is slow. When the mathematical relation 
of the motion of such moving objects can be identified, it is possible to estimate the possibility of a collision of such 
moving objects. 

(Background Art) 

50 

(Prior art) 

[0003] Some methods for identifying the mathematical relation between measured values are proposed. For exam- 
ple, a least square root method for determining the factors a and b in a linear relation y = ax + b is widely employed, 
55 under the assumption that a linear relation y = ax + b stands between the input data. Also a least square root method 
after a logarithmic transformation is widely employed. 

[0004] Methods for identifying the mathematical relation between the measured data are explained compactly in 
"Statistics for analytical chemistry" by J. C. Miller / J. N. Miller, which is translated into Japanese by Munemori Makoto 
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and published in Japan from Kyouritu Shuppan in 1991 . 

[0005] Japanese patent application JP-5-334431-A discloses an apparatus for giving a mathematical function, 
which approximates data of points on a line. This apparatus is not applicable, when the data is not data of points on a 
line, or when the data is data of over three dimension. 
5 [0006] Japanese patent application JP-5-266063-A discloses an apparatus for interpolating data in n dimension 
space, using a super surface in a (n+1 ) dimension space. This apparatus can interpolate data, but does not identify the 
mathematical relation of the data. 

(object of the Invention) 

w 

[0007] An object of the present invention is to propose an apparatus and method for outputting a mathematical rela- 
tion between base variables (x1 , x2, ...,xp), when a set of input data d is comprised of p base variables (x1 , x2, ...,xp), 
and a plurality of data sets d(i) of such a data set d are inputted, where the symbols V is a parameter for distinguishing 
the sets of data set. 

15 

( Disclosure of Invention ) 
(glossary) 

20 [0008] The meaning of symbols and terms used in this Specification and the Claims are explained, before explain- 
ing the present invention. 

Symbol " A " is a power (for example, (-2) A 3 = -8). Power is expressed, also using a suffix (for example, (-2) 3 = -8). 
Symbol M is a multiplication (for example 2*3 = 6). However, this symbol is abbreviated, when the meaning is 
25 obvious (for example, 2x + 3y = 0 means 2»x + 3»y = 0). 

Symbol " ( a • b ) " is an inner product (scalar product) of vectors. 

"Base variables x1 ,...,xp" are names of variables of data inputted into the apparatus according to the present inven- 
tion (for example coordinates x, y, z pressure p, or time t). 

- "Input data d" is a set of base variables (x1 , x2, ... xp). 

30 - "d = (x1 xp)" means that the number of the base variables of the input data d is p, and the final base variable is 

xp. 

- "Input data d(i)" is the i-th set of the input data d. Symbol V is called "data specifying parameter", 
"xji" is the value of the j-th variable of the i-th set of data d(i). 

"Mathematical function program" is a program for outputting a value corresponding to input reference values, after 
35 executing a mathematical calculation or using a table. 

- "An input reference and a function specifying reference": For example, a function exp(k • x) can be considered as a 
two variable function, in such a case, the function exp(k • x) has two input variables k and x. However, this function 
can be considered also as a one variable function exp(k • x). In this case, the function has an input reference x and 
a reference k for specifying the function. References for specifying the function are called "function specifying ref- 

40 erence", and references as input values are called "input reference". 

"Function specifying parameter m" is a parameter for specifying a mathematical function program gm. 
"Base function" is a function gk in a set of mathematical function programs gm stored in a mathematical function 
program storing memory, and specified by the function specifying parameter (m = k). The input variables to be 
inputted to the input reference of the mathematical function gk are specified. The function specifying references are 

45 specified when it is necessary. When the function is a constant function, which outputs a constant value irrespec- 
tively to the input, it is not necessary to specify the input references. 

- "Candidate mathematical relation" is a set of base functions (f1 , f2, fq). In the present invention, a mathematical 
relation is identified as a linear combination of base functions. 

"A temporal candidate mathematical relation" is a set of base functions selected temporarily as a candidate math- 
50 ematical relation, when the most appropriate mathematical relation is sought, by changing the candidate mathe- 
matical relation. 

"(f1 fq)" means that the candidate mathematical relation is comprised of q base functions, and the final base 

function is fq. 

"Vector component Fki" is a value of a base function fk, corresponding to the values of the base variable (x1 , .... xp) 
55 of the i-th input data d(i) inputted to the input reference of the base function fk. 

"Vector F(i)" is a vector in a q dimensional space comprised of vector components Fki (k = 1 q). 

- "Vector component array VE_ARRAY" is a two dimensional array storing vector components Fki. 

- "Vector space P is a space spanned by q dimensional vectors F(i). 
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- "Plane 11 is a set of points in a q dimensional vector space, which satisfies the following relation, where the LI , L2, 

Lq are constant (Appendix B): 

L1 -x1 + ... + Lq-xq=0 (1) 

5 

- "Dimension cosine" is a set of numerals (L1 , Lq) which satisfy the expression (1 ), and normalized as follows: 

(L1) 2 + (L2) 2 + ... + (Lq) 2 = 1 (2) 

w - "Mapping plane" is a mapping of a set of input data d(i) into the vector space F f and the mapping is a plane or nearly 
a plane. Also, a plane the square sum of the perpendicular lines from mapped points to the plane is zero or nearly 
equal to zero (namely, smaller than a predetermined valueXO), is called a mapping plane, (Appendix A). 

- "Plane degree index PI" is an index indicating in what degree the vectors F(1 ),..., F(q) can be regarded as a plane. 
Plane degree index can be defined on any industrial or commercial ground (Appendix D). 

15 - "Correlation sum < alb > " is a sum of products of a vector component Fai obtained from a base function fa and a vec- 
tor component Fbi obtained from a base function fb, the sum is carried out in respect with the input data specifying 
parameter i (Appendix F). 

felb>= £(Fai •Fbi) (3a) 

20 

The term (alb) means also the sum of the products multiplied by a square of weight Wi 2 : 

<alb)=ZWi 2 (Fai-Fbi) (3b) 
25 - "Correlation sum matrix C" is a matrix, the elements of which are correlation sums, 
(means to attain the object) 



[0009] The object of the present invention is attained by the apparatus for identifying the mathematical relation 
30 according to claim 1, and the method for identifying the mathematical relation according to claim 10. The principle of 
the present invention is explained in Appendix A. 

[0010] A series of data d(i) comprised of p basic variables xj (j = 1 , .... p) are acquired by a data acquisition means. 
The unit of data is d(i) = (x1i, x2i, xpi) . 

[0011] A plurality of mathematical function programs gm are stored in a mathematical program storing memory. 

35 [0012] A base function defining means defines a base function by specifying a specific mathematical function pro- 
gram gk in the mathematical function program gm in the mathematical program storing memory, by specifying the 
parameter (m = k ), and further specifying the base variables to be inputted into the input reference of the mathematical 
function. Also mathematical function specifying parameters are specified, when it is necessary. 
[0013] A candidate mathematical relation specifying means specifies a temporal mathematical relation by specify- 

40 ing a set of base functions (f 1 , fq). In this invention, the mathematical relation is identified as a combination of those 

base functions (f 1 fq). It is possible to design the candidate mathematical relation specifying means to specify a 

series of temporal candidate mathematical relations sequentially. 

[0014] A vector component acquisition means sends the values of the basic variables to oath of the input refer- 
ences of base functions fk (k = 1 f .... q) included in the candidate mathematical relation, for each data set d(i). The vec- 
45 tor component acquisition means gets the output value Fki of the base functions fk and stores them into a two 
dimensional array. The array is called VC_ARRAY. When the base functions are fk (k = 1 , .... q) t and s sets of data are 
sent, the vector component may VC_ARRAY is a two dimensional array of (q x s) size. 

[001 5] A set of Fji is considered as a vector F(i). And the values of the functions are stored in the vector component 
array VC_ARRAY for every vector F(i), as follows: 

50 

F(1) = (F11,F21 Fq1) 

F(2) = (F12, F22 Fq2) 



55 [001 6] When q sets of data in the array VE_ARRAY are considered as a coordinate of points P(i) in a q dimensional 
space F, and if the combination of the candidate mathematical relation is correct, the points shall be found in a plane 
(mapping plane) in the vector space F. If the candidate mathematical relation is approximately correct, those points shall 
distribute near to a plane (mapping plane). As explained in Appendix C, the mapping plane can be determined from a 
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determinant of Fji, or can be determined as a plane, the square sum of perpendicular lines from the points P(i) to the 
plane is the minimum. 

[0017] A direction cosine acquisition means calculates and outputs the direction cosine (L1 Lq) of a plane on 

which the points P(i) (i = 1 , .... q) are found, or a plane, the square sum of perpendicular lines from the points P(i) to the 
5 plane are the minimum. 

[0018] A mathematical relation outputting means identifies the mathematical relation as a combination of the base 
functions fk (k = 1 q), where the combination coefficients of the base functions are the direction cosine. The math- 
ematical relation outputting means outputs the mathematical relation in a form of the following linear combination or a 
mathematical relation deduced from it: 

10 

q 

£(Lk-fk) = 0 (4) 

k=1 

15 

[0019] According to an embodiment, the direction cosine acquisition means acquires the direction cosine, using a 
cofactor of a determinant V of a set of values of mathematical functions Fji. The round of this method is explained in 
Appendix C. 

[0020] According to another embodiment, the direction cosine acquisition means acquires the direction cosine from 
20 an eigenvector of a correlation sum matrix C obtained from an array of sets of the values of base functions. The ground 
of this method is explained in Appendix F. 

[0021 ] In another embodiment, the base function defining means defines a plurality of candidate mathematical rela- 
tions. And so far as the eigenvalue X calculated by the direction cosine acquisition means is not smaller than a prede- 
termined value, another candidate mathematical relation is set as a new temporal candidate mathematical relation to 
25 repeat to try to acquire a new eigenvalue, until the eigenvalue becomes smaller than the predetermined value. As a 
result, an appropriate mathematical relation can be acquired automatically. 

[0022] In another embodiment, the mathematical relation identifying apparatus is incorporated in an apparatus for 
interpolating or extrapolating data. When an empirical formula of a system includes a plurality of parameters is required, 
a candidate mathematical relation of a combination of a base functions is assumed, on an empirical or theoretical 
30 ground. The empirical formula can be identified, by acquiring the direction cosine, using the apparatus according to the 
present invention and estimating the validity of the expression. 

[0023] In another embodiment, the apparatus for identifying a mathematical relation is incorporated in an apparatus 
for measuring a trajectory. At first, infinite sets of the coordinates of points on a trajectory of a moving object (for exam- 
ple, golf club head) is acquired. Previously a candidate mathematical relation is assumed on an empirical or theoretical 
35 ground. Then the direction cosine is acquired, using the apparatus according to the present invention. An empirical for- 
mula of a trajectory can be obtained from the direction cosine, after estimating the validity. 

[0024] In another embodiment, the apparatus for identifying a mathematical relation is incorporated in an apparatus 
for measuring exterior or interior form of an object, (for example, face, foot, or shoe). At first, coordinates of infinite sets 
of points on the surface of an object is acquired. Previously a candidate mathematical relation is assumed on an empir- 
40 ical or theoretical ground. Then, the direction cosine is acquired, using the apparatus according to the present invention. 
An empirical formula representing the form can be obtained from the direction cosine, after estimating the validity. 
[0025] The present invention can be realized as a memory media, which can be read out by a computer, in which 
the main part of the invention according to claim 1 is coded as software. 

[0026] An embodiment of the present invention is a method corresponding to the apparatus for identifying a math- 
45 ematical relation according to claim 1. 

(Brief Description of Drawings) 

[0027] 

50 

Fig.1 is a block diagram of an apparatus for identifying a mathematical relation according to one of the best mode 
of the present invention. 

Fig. 2 is a flow chart showing the flow of the procedure in the apparatus and method for identifying a mathematical 
relation according to the present invention. 
55 Fig. 3 shows an application for identifying the form of a human foot. 

Fig .4 shows an application for identifying the form of the interior of a shoe. 
Fig.5 shows an application for identifying the trajectory of a golf club head. 
Fig.6 shows an application for identifying the trajectory of a ball thrown by a pitcher. 
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Fig.7 shows an application for identifying a relation between a pressure P, temperature T, and concentrations X, Y, 
and Z. 

(Best mode for Carrying out the Invention ) 

5 

[0028] Figs 1 and 2 show the apparatus and the method of the present invention. 

A. Data acquisition means 1 and data acquisition procedure P1 : 

w [0029] A data acquisition means 1 in the apparatus for identifying a mathematical relation acquires data from a (not 
shown) data outputting means in a preceding stage. Examples of the data outputting means are a keyboard of a com- 
puter for inputting a numerical data; a position coordinate specifying device for specifying a coordinate of a position 
specified by a cursor and/or mouse device in a display of a computer; a position data acquisition means for acquiring 
data from stylus for measuring a position coordinate of a point on the surface of a model object in a teaching-in type 

is working machine; a position coordinate detecting means, which gets image information of a moving body using a CCD 
camera and transforms the image information to data of a position coordinate; and a position coordinate detecting 
means for acquiring position coordinate of a moving body using a radar. 

[0030] The data acquisition means acquires a series of data d(i). Each data d(i) is comprised of p basic variables 
xj (j=1 , p). Figs. 3-7 show examples of data d(i) to be acquired by the data acquisition means. 

20 [0031] Fig. 3 (a), (b) and (c) show examples of foot. In this example, 2 dimensional coordinates of points P1 1 .... P1 4 
are measured for specifying the form of the foot. In this case, each data d is comprised of a two basic variables (x, y) 
(d= (x, y) ). If it is necessary, the data acquisition means can have a coordinate transforming means for transforming an 
orthonormal coordinate (x, y) to a polar coordinate (r, 9), the origin of which is a specified point, for example a point O 
in the figure separated distances Xq and y 0 from the most inner side and the most rear side of the foot, respectively. 

25 [0032] The data acquisition means (1) can be designed to measure a three dimensional coordinate (x, y, z) of a 
point of a foot. In this case, each data d(i) is comprised of three basic variables (x, y, z). If it is necessary, the data acqui- 
sition means (1 ) can have a coordinate transforming means for transforming the three dimensional coordinate to a cylin- 
drical coordinate (r, 9, z), the center line of which is a vertical line (z axis), for example. 

[0033] In the example shown in Fig. 4, two dimensional coordinates of points P1, P14 on an interior surface of 
30 shoe are measured. In this case, each data d(i) is comprised of a two basic variables (x, y) (d = (x, y) ). The data acqui- 
sition means (1) can be designed to measure a three dimensional coordinate (x, y, z) of points of a foot. In this case, 
each data d(i) is comprised of three basic variables. 

[0034] If it is necessary, the data acquisition means can have a coordinate transforming means for transforming an 
orthonormal coordinate (x, y) to a polar coordinate (r, 9), or to a cylindrical coordinate (r, 6, z), similarly as the case of 
35 Rg.3. 

[0035] In the example shown in Fig.5, coordinates of points (for example, 15 points) on a trajectory of a golf club 
head are measured as a function of time t. Each data d(i) is comprised of four basic variables (t, x, y, z). 
[0036] In the example shown in Fig.6, coordinates of points (for example, 5 points) on a trajectory of a ball thrown 
by a pitcher are measured as a function of time t. Each data d(i) is comprised of four basic variables (t, x, y, z). 

40 [0037] The example shown in fig.7 is a case to study a relation between a pressure P, temperature T, concentration 
X of a component A, concentration Y of a component B and concentration Z of a component C. Data (P, T, X, Y, Z) are 
measured at several fixed values of X and T. In the figure, some curved surfaces representing relations between P, X, 
Z at fixed values of the parameters X and T are shown. The data acquisition means 1 acquires data (P, T, X, Y, Z) at 
points d(1) d(20) on the curved surface. In this case, each data d(i) is comprised of five basic variables. 

45 [0038] The data acquisition means (1 ) can be realized as a file (F1 ) of a data acquisition program. 

B. Input data array forming means 2 and data array forming procedure P2: 

[0039] It is preferable to store the series of data d(i) in a data array DATA_ARRAY in an input data memory. The 

so data unit in it is a set of data d(i). When s sets of data d(i) (i = 1 s) are acquired, the data array DATA_ARRAY is a 

two dimensional array, and the size of the array is (p x x). For Figs. 3 and 4, s = 14. And for Fig.5, 6, 7, respectively, s= 
15, s=5ands = 20. 

[0040] The data array forming means can be realized as a file (F2) of a data array forming program. However, the 
forming of data array DATA_ARRAY is not inevitable. It is possible to design to send the data acquire by the data acqui- 
55 sition means 1 directly to a vector component acquisition means 6. which will be explained hereinafter. 
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C. Mathematical function program forming means 3 and mathematical program storing procedure P3: 

[0041 ] Mathematical function program are formed by combining a general purpose mathematical function and tem- 
poral input references IN(1), IN(2), ...etc. The general purpose mathematical functions include, for example functions 

5 for, addition, subtraction, multiplication, division, power, exponential function, logarithm function, trigonometric function, 
and a specially programmed functions (for example, BessePs function). Depending on the mathematical function, some 
mathematical function specifying reference shall be specified for determining a mathematical function program. 
[0042] The mathematical function programs will be stored in a mathematical function program memory. The 
number of the temporal input reference variables depends upon the mathematical function. 

10 [0043] The mathematical function program can be realized as a file F3 of list of formed program list or a program 
for forming a mathematical program. 

[0044] For example, a mathematical function program g = IN(1 ) • IN(2) can be defined, by combining a multiplica- 
tion function incorporated in a computer and two temporal input reference variables IN(1) and IN(2). 
[0045] Some other examples are shown below: 
15 Constant function has no input reference. 

gO = 1 (5-0) 
[0046] Examples of mathematical function having only one input reference variable (IN(1)): 

20 

gi = iN(i) (5-1) 

g2 = IN(1) 2 (5-2) 

25 g3 = IN(1) 3 (5-3) 

g4 = cos(IN(1)) (5-4) 

g5 = cos(2-IN(1)) (5-5) 

30 

g6 = sin(IN(1)) (5-6) 

g7 = sin(2-IN(1)) (5-7) 

35 g8 = IN(1) •cos(INd)) (5-8) 

g9 = IN(1) -sin(IN(1)) (5-9) 
[0047] Examples of mathematical function having two input reference variables (IN(1) and IN(2)): 

40 

g10 = IN(1) -IN(2) (5-10) 

g11 =IN(1) -cos(IN(2)) (5-11) 

45 g12 = IN(1) -sin(IN(2)) (5-12) 

g13 = IN(1) 2 -IN(2) (5-13) 
[0048] Examples of mathematical function having three input reference variables (IN(1), IN (2) and IN(3)): 

50 

g14 = IN(1) *IN(2) -^(3) (5-14) 
[0049] Examples of mathematical function having four input reference variables (IN(1), IN (2) and IN(3)): 

55 g15 = (IN(1)-IN(2))/(IN(3)- IN(4)) (5-15) 

g16 = (IN(1)*IN(2)) • (IN(3) A IN(4)) (5-16) 
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[0050] In the explanation hereinafter, the functions gO, g1 , g2, g3, ... are used in a meaning as defined here. How- 
ever these definitions are no more than examples for explanation of concrete examples. 

D. Base function defining means 4 and base function defining procedure P4: 

5 

[0051] A base function defining means 4 defines a base function f, by specifying a function specifying parameter 
(m = k ), which specifies a mathematical function gk from functions gm in the mathematical function program memory 
(3), and further specifying a basic variable to be inputted into the r-th temporal input reference variable IN(r) of the spec- 
ified mathematical function gk. When the specified mathematical functions is a constant function g0, it is not necessary 
w do specify the input reference variable. Even when identical mathematical function is used, when basic variables to be 
inputted to its input reference variables IN(r) are different, they form different base functions. 

[0052] The base functions are stored in a base function defining may BF_ARRAY. When a set of input data d is 
comprised of basic variables (x1 , x2, x3), base functions corresponding to, for example, a constant function, function 
x1, function x2, function x1 *x2, function x1 *x3, function (x1)2, function (x2)2, and differential dx1/dx3 are expressed 
15 as follows, using the mathematical function programs gO, g1 , 



BF_ARRAY j 


output 


mathematical function 
program 


basic variables to be 
inputted to temporal 
input references 


1 


" go 






X1 


gi 


IN(1): 


= x1 


X1 2 


g2 


IN(1) = 


= X1 


x1 »x2 


g8 


IN(1) = 


= x1 






ln(2) = 


= x2 


x1 *x3 


ge 


IN(1): 


= x1 






ln(2) = 


= x3 


dx1/dx3 (differential) 


gi4 


IN(1): 


= x1(i+1) 






IN(2) = 


= x1(i) 






IN(3): 


= x3(i+1) 






!N(4) : 


= x3(i) 



40 

[0053] The base function defining means 4 can be realized as a file F4 storing a program for defining base functions 
or a memory storing defined base functions. 

E. Candidate mathematical function specifying means 5 and candidate mathematical functions specifying procedure 
45 P5: 

[0054] A candidate mathematical function specifying means 5 specifies a set of base functions selected from the 
base functions defined by the base function defining means. The temporal candidate mathematical relation shall be 
specified according to a theoretical ground or an empirical law. The temporal candidate mathematical relations can not 
so be determined uniquely, but shall be determined by a person skilled in the art so that the mathematical relation of the 
object can be appropriately expressed, using a small number of simple mathematical functions. 
[0055] For example, following mathematical relations can be supposed for representing a trajectory of the golf club 
head shown in Fig. 5: 

55 #1)p*x 2 +q*y 2 + r- z 2 + a*x + b»y + cz + dO+d1't + d2M 2 =0 (6) 

#2)a*x + b*y + c- z + d»1=0 (7) 
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[0056] The expression #1 ) represents a model of a movement on an ellipsoid, which has a fixed center and the size 
is changing according to a quadric function of time t. The expression #2) represents a model of a movement on a plane 
in an actual space, the distance from the origin to the plane is d, and the direction cosine of the plane is (a, b, c). When 
it is estimated that both of the mathematical relations stand, the trajectory of the golf club head can be identified as a 
5 cross sectional line of a size changing ellipsoid and a static plane. 

[0057] In the case #1), the candidate mathematical relation can be defined by selecting the following set of base 
functions (f1 f9): 



CF(1) 


Base function 


Mathematical function 
program 


Reference Variable IN(1) 


f1=x 2 


92 


X 


f2 = y 2 


92 


y 


f3 = z 2 


g2 


z 


f4 = x 


gi 


X 


f5 = y 


gi 


y 


f6 = z 


gi 


z 


f7 = 1 


go 




f8 = t 


gi 


t 


f9 = t 2 


g2 


t 



[0058] In the case #2), the candidate mathematical relation can be defined by selecting the following set of base 
30 functions (M, .... f4): 



CF(2) 


Base function 


Mathematical function 
program 


Reference Variable IN(1) 


f1 =x 


gi 


X 


f2 = y 


gi 


y 


f3 = z 


gi 


z 


f 4 = 1 


go 





45 [0059] When input data d(i) of the form of a foot shown in Fig.3 are expressed in terms of a polar coordinate (r, 0), 
in reference with an origin O, and when the form of the foot is approximated as forms (d), (e), and (f) shown in Fig.3 t 
namely, the forms d, e and f are represented in a form of r = r(6) , as shown in (g) of Fig.3, an example of the candidate 
mathematical relation between r and 9 is as follows: 

so #1) r = a + b • 9 + c • 0 2 + d • cose +e • cos20 + f • sine + g • sin2e + h • Ocose + i • 8sine (8) 

[0060] In this case, the candidate mathematical relation can be defined by selecting the following set of base func- 
tions (M f10): 
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5 



10 



20 



CF(3) 


Base function 


Mathematical function 
program 


Reference Variable IN(1) 


f1 = 1 


go 




f2 = e 




e 


fa a2 
to = D 


9 2 




f4 = cosG 


94 


e 


f5 = cos2G 


95 


e 


f6 = sine 


96 


e 


f7 = sin2e 


97 


e 


f8 = ecose 


98 


e 


f9 = esine 


99 


e 


f10 = r 


911 


r 



[0061] This set of base functions are no more than an example. It is possible to select another set of base functions, 
25 for specifying another candidate mathematical relation. In connection with the form of the interior surface of a shoe 
shown in Fig.4. It is possible to specify a candidate mathematical relation, by selecting a similar set of base functions. 
[0062] In connection with a trajectory of a ball thrown by a pitcher shown in Fig.6 t a candidate mathematical relation 
in an orthonormal coordinate system (x, y, z) can be assumed as follows: 

30 x + a0 + a1 •t + a2*t 2 = 0 (9-1) 

y + b0 + b1 -t + b2*t 2 = 0 (9-1) 

z + cO + d M + c2-t 2 = 0 (9-1) 

35 

[0063] In this case, the candidate mathematical relation between x and t can be defined by selecting the following 
set of base functions (f1 , f4): 



50 



CF (4) 


Base function 


Mathematical function 
program 


Reference Variable 1N(1) 


H =x 


91 


X 


f 2 = 1 


go 




f3=t 


gi 


t 


f4 = t 


g2 


t 



[0064] Also in connection with the y, z, similar temporal candidate mathematical relations can be specified. 
[0065] In connection with Fig.7, assuming that there is a following empirical relation between the pressure P, tem- 
perature T, concentration of a composition A, concentration Y of a component B, and a concentration Z of a product: 

55 

Z/XY=K0 + K1 •P + K2-T+K3-X (10-1) 
[0066] This expression can be rewritten as follows: 
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K0 • XY + K1 • XYP + K2 • XYT + K3 • X 2 Y - Z = 0 (10-2) 

[0067] In this case, a candidate mathematical relation can be defined by selecting the following set of base func- 
tions (f1. f 5): 

f1 = X-Y, f2 = X-T-P, f3 = X-Y-T, f4 = X 2 -T, f5 = Z 



CF (5) 


Base function 


Mathematical function 


Reference Variable 




program 












IN(1) 


IN(2) 


IN(3) 


f1 =XY 


gio 


X 


Y 




f2 = XYP 


g14 


X 


Y 


P 


f3 = XYT 


gi4 


X 


Y 


T 


f4 = X 2 T 


gi3 


X 


T 




f5 = Z 


gi 


Z 







[0068] A candidate mathematical relation specifying means can be realized as a file F5 storing a program for spec- 
ifying a candidate mathematical relation and/or as a memory storing sets of base functions as candidate mathematical 
relations. This program can include a program for sequentially specifying different candidate mathematical relations 
according to a signal from a validity estimating means. 

F. Vector component acquisition means 6, vector component acquisition procedure P6 and memory 7 for vector com- 
ponent array: 

[0069] Value xji of basic variable xj is sent to the input reference variable IN(r) of each of the base functions, so as 
to be substituted into the mathematical program, then the acquired value Fki of the function is outputted as a vector 
component. 

[0070] It is preferable that the vector component acquisition means (6) comprises an array INA for temporary stor- 
ing the values of basic variables transferred from a data array DATA_ARRAY in the memory (2) for input data, and that 
the sets of data d(i) are sent form the array DATA_ARRY for every input data specifying parameter i to this array. 
[0071] It is further preferable that the values xmi (i = 1 , 2, 3,...) of the m-th base function are sequentially transferred 
into the m-th element INA[m] of this array INA. For example, when d is a set of data X, Y, Z, P t T, namely d = (X, Y, Z, P, 
T), data of X, Y ( Z, P, T are stored in the element INA[1], INA[2], INA[3], INA[4], and INA[5] of the array INA. It is prefer- 
able that the base function defining means defines the content of the input reference variable of a base function, by 
defining what variable shall be stored in each of the element INA[m] of this array INA. 

[0072] For example, the input reference variables of the base functions f1 (= XY), f2 ( = XYP), f3 ( = XYT), f4 ( = 
X 2 ^ and f5 ( = Z) in CF(5) are variables X, Y, Z, P and T. And the base functions can be defined as follows: 



CF(6) 




Mathematical function 


Input reference variables 




program 












IN(1) 


IN(2) 


IN(3) 


f1 


g10 


INA[1] 


INA[2] 




f2 


gl4 


INA[1] 


INA[2] 


INA[4] 


f3 


gi4 


INA[1] 


INA[2] 


INA[5] 
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(continued) 



CF(6) 




Mathematical function 


Input reference variables 




program 












IN(1) 


IN(2) 


IN(3) 


f4 


g13 


INA(1] 


INA[5] 




f5 


gi 


INA[3] 







[0073] When data d(i) are transferred to the input variable array INA for every input data specifying parameter i, 
basic variables xr are transferred to the corresponding input reference variables IN(r) of each of the base functions 
through this array INA. 

[0074] The values xji of the basic variable xj are transferred to the input reference variable of each of the base func- 
tion, indirectly through the array INA or directly without passing through this array. Then the values are substituted into 
the mathematical function program and the values Fki are acquired. The values Fki are stored into a vector component 
may VC_ARRAY. By the way, it is not necessary to acquire the values Fki for all of the data set d(i). An embodiment to 
acquire the values Fki for only a part of the data sets d(i) is possible. 

[0075] When a set of data is comprised of basic variables X, Y, Z, P and T, namely d = (X, Y, Z, P, T), data d(i) f for 
example, d(5) is transferred to the input reference variabls IN(r) of the base functions f 1 = XY , f2 = XYP , f3 =XYT , 
f4 s= X 3 T , and f5 = Z of temporal candidate mathematical relation #1) through the array INAfr] as follows: 



INA[1] 


= X5-> IN(1) 


= X5 


INA[2] 


= Y5 -> IN(2) 


= Y5 


INA[3] 


= Z5 -» IN(3) 


= Z5 


INA[4] 


= P5-> IN(4) 


= P5 


INA(5] 


= T5 -> IN(5) 


= T5 



[0076] Then the values of each of the base functions are acquired: 

F15 = XY = g10(X5, Y5) = X5 • Y5 
F25 = XYP = g1 4(X5, Y5, P5) = X5 • Y5 • P5 
F35 = XYT = g1 4(X5, Y5, T5) = X5 • Y5 • Y5 
F45 = X 2 Y = g13(X5, Y5) = (X5) 2 • Y5 
F55 = Z = g1 (Z) = Z5 

[0077] Vector components (F1i, F2i, F3i, F4i, F5i) for each data specifying parameter i will be acquired similarly. 
The vector component acquiring means 6 can be realized as a file F6 storing a vector component acquisition program 
for carrying out the above mentioned functions. 

[0078] A set of the function values (F1i, F2i, .... Fqi) corresponding to a specific value of the parameter i forms a 
vector F(i) in a q dimensional space. The vectors F(i) (i = 1 , 2, ...) are stored in a vector component may VC_ARRAY in 
a memory for a vector component array. The array is a two dimensional array having q columns. An array VC_ARRAY 
is shown below, where s is the number of data sets: 



VC_ARRAY 


F11 


F21 


F32 




Fq1 
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(continued) 



5 



VC_ARRAY 


F12 


F22 


F32 




Fq2 


F1s 


F2s 


F3s 




Fqs 



G. Direction cosine acquisition means 8 and direction cosine acquisition procedure P8: 

w 

[0079] When the vectors F(i) stored in the vector component array VC ARRAY form a plane in a q dimensional vec- 
tor space F, or they distribute near to a plane in the space, the direction cosine acquisition means (8) acquire the direc- 
tion cosine of the plane, using a calculation of cofactors of a determinant or by obtaining an eignevalue of a matrix, for 
example. 

15 

G1. Direction cosine acquiring means using cofactors of a determinant: 



20 



25 



30 



[0080J The base of the mathematical ground this method is explained in Appendix C. This method is suitable in 
case that the number of the data sets is small. 

[0081] The direction cosine acquisition means using a cofactor of a determinant in a best mode comprises: a deter- 
minant array forming means 8A1 , a cofactor calculating means 8A2, a plane degree index calculating means 8A3, plane 
degree judging means 8A4 and a direction cosine outputting means 8A5. 

[0082] A determinant array forming means 8A1 selects q vectors arbitrary from the array VC_ARRAY, and forms a 
determinant of (qxq). When q vectors F(1), .... F(q) are selected, the determinant is as follows: 



Fll 
F12 



F21 
F22 



F31 
F32 



Flq F2q F3q 



Fql 
Fq2 

Fqq 



-.(11) 



35 [0083] The cofactor calculating means 8A2 calculates minors As1 , As2, .... Asq for each row s of this determinant. 
[0084] The direction cosine outputting means output multiplies (-1 ) A (k+s) to each of the minors Ask, to obtain cofac- 
tors Bsk, and output the cofactors as direction cosines Lk. In this specification, the symbol " An means a power. When 
the square sum of the L1 , L2 Lq are normalized to be one, a direction cosine in an ordinal meaning can be obtained. 

40 (L1) 2 + (L2) 2 + ... + (Lq) 2 = 1 (12) 

[0085] The plane degree index calculating means 8A3 means calculates a plane degree index PI form the minors. 
The plane degree index is explained in Appendix D. 

[0086] The plane degree judging means 8A4 judges whether the difference between the plane degree index 
45 obtained by the plane degree index calculating means and a predetermined standard value pi is within a allowable 
region dPI or not. The standard value pi of the plane degree index PI and the allowable region dPI are explained in 
Appendix D. 

[0087] When the difference is out of the allowable region, the sets of the data may include sets of data having a 
large error, thus this sets of data shall be abandoned, and sends a signal to the determinant array forming means to 
so select other sets of data. When the difference is within the allowable region, the plane degree judging means sends a 
signal to the direction cosine outputting means (8A5) to output the cofactors or to output mean values of each of the 
cofactors (L1 >, .... <Lq>. 

[0088] As explained in Appendix E, the plane degree judging means (8A4) can have a function to judge whether 
specific cofactors (or minors) are always equal to zero or not. Another candidate mathematical relation shall be set, 
55 when some specific cofactors are always equal to zero. 

[0089] In a best mode, the apparatus according to the present invention calculates a plane degree index and esti- 
mate it, as has been explained. However, when it is clear that the errors of data d(i) are small, the calculation is not nec- 
essary. The direction cosine can be acquired, by calculating a set of cofactors for an arbitrary row of the determinant, 
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without passing through the procedure of calculating the plane degree index and judging the plane degree. The deter- 
minant correspond to vectors F(1 ), F(2), .... F(q). In this case, the normalization is not inevitable. 
[0090] The method for calculation of the determinant, minors, and cofactors using a computer are well known, thus 
the explanations about them are omitted. 

G2. Direction cosine acquiring means using eignevalue of a matrix: 

[0091] This means and method are suitable, when a lot of data sets d(i) can be acquired. The ground of this method 
is explained in Appendix R As explained there, an eigenvalue corresponds to the square sum the perpendicular lines 
from each of the vectors F(i) to a mapping plane. It is explained in Appendix G that all the eigenvalues are real number 
and non-negative. By obtaining a positive eigenvalue smaller than a value X0, which is a predetermined value of the 
square sum of the length of the perpendicular line, and obtaining the eigenvector corresponding to the eigenvalue, the 
mapping plane can be determined. 

[0092] In the best mode of the direction cosine acquisition means 8 using an eigenvalue comprises: a correlation 
sum calculating means 8B1 , a eigenvalue calculating means 8B2 and a eigenvector calculating means 8B2. 

[0093] The correlation sum calculating means 8B1 calculates a sum of the product of a-th (a = 1 q) element Fai 

and b-th (b = 1 q) element Fbi in the vector component array VE_ARRAY (Appendix (F3), (F17)). The obtained cor- 
relation sums form a correlation matrix C and are stored in a correlation sum array CO_ARRAY Namely, the element 
<alb) of the correlation sum is stored at a-th row, b-th column in the array CO_ARRAY. An example of the correlation 
sun array is shown below. 



[0094] A secular equation corresponding to the correlation sum C is formed. And the eigenvalue calculating means 
8B2 solves the secular equation, to obtain positive eigenvalues X1 , X2, ... of the secular equation, which are close to 
zero, namely smaller than a predetermined value X0. 

[0095] It is possible to obtain all the eigenvalues, and to select eigenvalues smaller than a predetermined value XO. 
In this case, all the eigenvalues must be obtained by an all eigenvalue calculation means ALL, and smaller eigenval- 
ues are selected using a comparator COMR However, when the secular equation is developed into a polynomial of (- 
X) using a polynomial forming means POLYNOM, the eigenvalues can be obtained easily in short time. In the develop- 
ment of the equation, t-th traces of the secular equation, which can be obtained using a trace calculation means, can 
be used. The zero points of the polynomial in a region [0, X0] are calculated. This method is explained in Appendix H. 
[0096] In the next step, the eigenvector calculating means 8B3 calculates an eigenvector corresponding to eigen- 
values Xm. The method to obtain an eigenvector from an eigenvalue is known, thus the explanation is omitted. 
[0097] By the way, the number of eigenvalues close to zero (eigenvalues smaller than a predetermined value XO) is 
not always one. For example, assuming that points P(i) in a q dimensional space correspond to data d(i). When the 
points P(i) distribute along a line, which is a cross sectional line of two (q-1) dimensional planes in the q dimensional 
space, one can understand that the data d(i) are mapped to two mapping planes. For example, assuming that a trajec- 
tory of a golf club head is a cross sectional line of an ellipsoid and a plane, when base functions corresponding to the 
expression (6) is selected as a candidate mathematical relation, two sets of eigenvectors, one of them corresponds to 
an equation of an ellipsoid and the other corresponds to an equation of a plane, can be obtained. 
[0098] The direction cosine acquisition means can be realized as a file F8 storing a direction cosine acquisition pro- 
gram, which acquires a direction cosine according to G1 or G2. 

H Validity estimating means 9 and validity estimating procedure P9: 

[0099] A validity estimating means 9 estimates the validity of the obtained mathematical relation on the basis of the 
eigenvalue, which corresponds to the obtained direction cosine acquired in the direction cosine acquisition means P9; 
[0100] For example, zero points in region [0, XO] can be obtained after developing the secular equation into a poly- 




Correlation sum matrix C 
<lll> <1|2> ... <1| 
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nomial. However, the allowance of the measurement error shall be determined independently, on the basis of industrial 
or commercial grounds. When a smaller allowance of the measurement error is permissible, the allowance of eigne- 
value Xrei shall be set smaller than the allowance value X0, so as to abandon the eigenvalues larger than the value A/ef. 
When the candidate mathematical relation is not be appropriate, or the sets of data may include data set including large 

5 error data, the eigenvalue will become larger than the value Xref. In such a case, the process returns to the procedure 
P6 to select another sets of data d(i) f or returns to the procedure P5 to select another candidate mathematical relation. 
[0101] Namely, when the eignevalue is not smaller than the predetermined value, the candidate mathematical rela- 
tion will be abandoned, and the process returns to the candidate mathematical relation specifying procedure P5, as 
shown in Figs.1 and 2. The candidate mathematical relation specifying means 5 selects another mathematical relation, 

w the vector component acquisition means 6 calculates the vector components, and the direction cosine acquisition 
means 8 acquires again a direction cosine. These procedures are repeated, until an appropriate eigenvalue is obtained. 
[01 02] When no candidate mathematical relation can be predicted, different temporal candidate mathematical rela- 
tions are set sequentially and/or different sets of data d(i) are set sequentially, and the mathematical relation identifica- 
tion procedure is repeated, so as to acquire an eigenvalue as small as possible. 

15 [0103] When a certain candidate mathematical relation can be predicted, for example, it is certain that the form is 
a paraboloid and only the identification of the values of the parameters is required, the base functions shall be fixed, so 
that the mathematical relation can be identified only by obtaining the direction cosine. It is no need to return to the can- 
didate mathematical relation specifying procedure P5. 

[0104] By the way, the validity estimating means is not inevitable to carry out the present invention. For example, 
20 when allowed time for data processing is short, and it is known that the error in the input data is small, the validity esti- 
mating means can be omitted. 

[0105] The validity estimating means 9 can be realized as a file F9 storing a validity estimating program for estimat- 
ing the validity. 

25 I. Mathematical relation outputting means 10 and mathematical relation outputting procedure P10: 

[0106] When the direction cosine of the mapping plane for each of the base function is Lk (k = 1 , 2, ...,q) ? and the 
validity estimating means 9 estimates that the eigenvalue is valid, mathematical relation outputting means 10 outputs 
an expression corresponding to the following mathematical relation as the acquired mathematical relation: 

30 

X( Lk - fk ) = ° ( 14 ) 

k=i 

[0107] For example, when the obtained direction cosine of the temporal candidate mathematical relation of the 
example of CF(5) is (L1 , L2, L3, L4, L5), the following mathematical expression or an expression or value deduced from 
the expression is outputted: 

40 L1 *f1 +L2-f2 + L3-f3 + L4*f4 + L5-f5==0 (15) 

[0108] Namely the following expression will be outputted: 



35 



45 



L1 • XY + L2 • XYP + L3 ♦ XYT + L4 • X 2 Y +L5 • Z = 0 (16) 

[0109] Because f1 = XY , f2 = XYP , f3 = XYT , f4 = X 2 Y , and f5 = Z , according to CF(5). 
[0110] The mathematical relation can be sent to a controlling apparatus, display, or printer, etc, which are con- 
nected to the apparatus of the present invention. The controlling apparatus will carry out a necessary calculation on the 
basis of the mathematical relation and sends a controlling signal to another apparatus. 
so [0111] The mathematical relation outputting means can be realized as a file F10 storing a mathematical relation 
outputting program for outputting a mathematical relation. 

[01 12] As explained, an apparatus for identifying a mathematical relation in best mode can be realized using: a file 
F1 storing a data acquisition program, a file F2 storing a data array forming program when it is necessary, a file F3 stor- 
ing prepared mathematical function programs and/or a mathematical program forming program, a memory storing base 
55 functions and/or a file F4 storing a base function defining program, a file F5 storing a program for specifying a candidate 
mathematical relation and/or specifying a memory storing sets of base functions as candidate mathematical relation, a 
file F6 storing a vector component acquisition program, a file F8 storing a direction cosine acquisition program, a file F9 
storing a validity estimating program, a file F10 storing a mathematical relation outputting program, and a file F1 1 stor- 
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ing a control program for controlling these programs to function as described in this specification. 
< Industrial Applicability ) 
5 [0113] 

A. The present invention can be applied as an apparatus for measuring a form of an object: For example, a form of 
a foot or a form of an interior of a shoe can be identified using the apparatus according to claim 1 , then the form of 
such an object can be specified using a set of direction cosine or a set of values deduced from the direction cosine. 

w Namely, selecting an appropriate set of base functions, a form can be expressed as a linear combination of the 
base functions. The form can be represented by the combination coefficient. A data base of a form can be con- 
structed, using this coefficient as a parameter. 

B. The present invention can be used as an apparatus for identifying a trajectory of a moving object. The mathe- 
15 matical relation of a trajectory of a golf club head, or a ball thrown by a pitcher, etc., can be identified, so that the 

trajectory can be analyzed or the trajectory in the future can be predicted. An application to an identification of a 
trajectory of a golf club head and its analysis are explained, in Appendix I. 

C. The present invention can be applicable as an analyzer of measured data. A mathematical relation can be iden- 
20 tified with relatively small number of experiments, even when a set of data is comprised of a lot of parameters, for 

example, a data of a chemical process. 

[01 14] There is a case that a differential equation can be identified from data, when the data shall be described by 
a differential equation. An example is shown in Appendix J. 

25 
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<Appendix> 

A. <Relation between combination coefficients and direction cosine> 

In the present invention, a mathematical relation is expressed as a linear 

combination of base functions fk (k = 1. 2, .... q). 

The principle of the present invention is that a linear combination of base 
functions is supposed as a candidate mathematical relation, and input data 
are transformed to sets of data in a q dimensional space spanned by the 
base functions, each set of data forms a vector in the space, when the vectors 
are found in a plane in the q dimensional space, the supposed mathematical 
relation is understood to be correct The essential idea of the present 
invention is to propose an apparatus and method for finding the direction 
cosine of the plane and for outputting the mathematical relation. 

For example, when sets of input data with two variables d(i) = (x(i), y(i)) (i = 
1...) are given, and the following mathematical relation is supposed to stand: 

a-x 1 + b-x + c-y + d-1 = 0 ... (Al) 

The sets of data (x, y) shall form a parabola curve. Although, when the data 
(x, y) are transformed by a set of functions, xl = x*, x2 = x, x3 = y, x4= w 
where (w =1), the following relation will be obtained: 

a-xl + b-x2 + coc3 + d-x4 = 0 ... (A2) 

x4=l ,..(A3) 

We introduce a 4 dimensional vector P(xl, x2, x3. x4) having vector 
components (xl, x2, x3, x4). A point P(i) corresponding to a set of data d(t) 
can be found at a cross sectional line of two planes, a first plane 1 satisfies 
the relation (Al), and a second plane 2 satisfies the relation (A3). 

The proportion of the direction cosine of the plane 1 are (a : b : c : d). 
Therefore, if one can determine the direction cosine of the plane 1 from the 
sets of input data d(i), and when the relation (A3) is combined, one can 
identify the mathematical expression to be a parabola curve. 

The present invention generalizes this idea. Namely at first, a set of base 
functions, inputs variables of which arc measured values, are defined. A 
candidate mathematical relation is predetermined by selecting a set of the 
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base functions. The measured valued are inputted into each base function to 
obtain the value of the base function. The set of the values of the base 
functions forms a vector space, namely the measured data are mapped into a 
vector space. When the mapping is a plane in the vector space, the direction 
cosine of the plane is calculated. Finally, the mathematical relation is 
identified as a linear combination of the base function. The combination 
coefficients of the linear combination are proportional to the direction cosine. 

When the mapping form a plane or substantially a plane in the vector space, 
in this description of the invention, we call it "a mapping plane". The term 
"mapping plane" is used also as a plane which gives the value zero or the 
minimum value of the sum of squares of distances between the plane and the 
mapped points of the data. 

When the mapped points do not form a mapping plane or cannot be 
approximated as a mapping plane, the data or the candidate mathematical 
relation is judged to be improper. Another set of data or another candidate 
mathematical relation will be selected, until a mapping plane can be found 
so as to obtain a proper mathematical relation between the data. 

B. <A plane in a vector space at higher dimension> 

In this specification, a "plane" is a set of points in a q dimensional vector 
space (xl, x2, .... xq), the points in the set satisfy the following relation: 

Ll-xl + ... + Lq-xq = 0 ... (Bl) 

where the Ll, L2, ... Lq are constant coefficients. 

Also a set of points on a intersecting line of a plurality of planes in a q 
dimensional vector space is called a "space". For example, a set of points on 
the crossing of two planes in a three dimensional space is a line, however, in 
this specification, it is called a "space". 

For example, when coordinates of an object (x, y) is represented by the 
following expression, which are functions of time t: 

x(t) = clt* + dl-t + cl ...(B2) 

y(t) = c2-t* + d2t + e2 ... (B3) 

When base functions xl to x5 defined as xl = x. x2 = y. x3 = t 2 , x4 = t, x5 = w, 
(w=l), are introduced, these relations can be expressed by a set of three 
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planes in a vector space (xl, x2 v x3, x4, x5): 
a-xl + clx3 + dl-x4 + el-x5 = 0 
a-x2 + c2-x3 + d2ot4 + e2-x5 = 0 
x5 = w=l 

On the other hand, a form of a surface can be expressed using a single 
mathematical expression. A form of a surface is not limited to a form of an 
object, for example, an outer surface of a foot, or an inner surface of a shoe, 
but also includes a surface in an abstract space, which expresses a 
mathematical relation between a plurality of parameters, for example, a 
relation between temperature T, pressure P, volume V and concentration C. 

C. <Relation between a cofactor and a direction cosine> 
When q vectors F(s) = (Fls, F2s, .... Fqs) (s = 1, 2, .... q) in a q dimensional 
space and the origin (0, 0) are included in a plane of a q dimensional 
space, the following determinant V is equal to zero: 
Fll F21 ... Fql 



F12 F22 ... Fq2 
... ... ••• 

Flq F2q ... Fqq 



(CI) 



For example, when base functions are fl, f2, 13, and the value of the base 
functions fl, f2, f3 corresponding to the detected values d(l), d(2), d(3) are 
(Fll, F21, F31), (F12, F22, F32), (F13, F23, F33), the points PI, P2, P3 
corresponding to these vectors and the origin (0,0,0) form a parallelogram. 
However, when the parallelogram is contained in a plane, the following 
determinant V is equal to zero. 
Fll F21 F31 

= 0 



V = 



F12 F22 F32 
F13 F23 F33 



It is known that the equation corresponding to the plane of the mathematical 
expression CI can be expressed by a determinant, in which one vector in the 
vectors F(l), F(2), ... F(q) are replaced to coordinate variables (xl, ... xq). For 
example, by replacing the q-th vector F(q) to (xl, x2.. . . ,xq). an equation of the 
plane can be obtained, as follows: 
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0 = 



Fll 
F12 



F21 
F22 



Fql 
Fq2 



Flz F2z ... Fqz 
xl x2 ... xq 
where (q-1) is replaced to z, namely (z 
understanding of the determinant. 



... (C2) 



= q-1), for facilitating the 



When this determinant is expanded in respect to the q-th row, the following 

linear combination will be obtained (Aqk is a minor, and s(q, k) is a sign), 
0 = s(q,l) Aql-xl + s(q. 2) Aq2-x2 + ... + s(q, q) -Aqq xq ... (C3) 
s(q,k) = (-l) A (q + k) ...(C4) 

This expansion can be rewritten as follows: 
0 = Bql -xl + Bq2x2 + ... + Bqq -xq ... (C5) 

where Bqk is a cofactor ( s(q, k) Aqk). 

On the other hand, direction cosine Ll Lq of a plane satisfies the 

following relation: 

0 = Ll-xl + ... 4- Lq-xq ... (C6) 

Consequently the following proportional relation stands: 
Ll: L2: ... Lq 
= Bql: Bq2: Bqq ... (C7) 

Therefore, the direction cosines is proportional to cofactors. 

When the determinant V = 0, the mapping is a plane, and the obtained 
proportion coefficients of the direction cosines is constant, irrespective to the 
row F(j) 0 = l f .». q), which is replaced in calculation of minors. 

When V^O, it is judged whether the vectors are found nearly in a plane or 
not, on the basis of an index PI. The index PI indicates in what degree the 
vectors F(l), ...» F(q) can be regarded as a plane. It is judged whether the 
index is within a predetermined region dPI or not. In case that V^O and 
that the index PI is found in a predetermined region dPI, it can be judged 
that the mapping is nearly a plane, namely, a mapping plane. 

The index PI and the predetermined region will be explained in Appendix D, 
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D. <Plane degree index PI and allowance dPI> 

A vector F(k) in a q dimensional vector space can be expressed as (Flk, 
FZk, fqk). When the origin (0. 0, 0) of the vector space and q vectors 

F(l) F(q) are in a common plane, the determinant IF! corresponding to 

the vectors F(l), F(q) is zero; 



IFI = 



Fll F21 ... Fql 
Flq F2q ... Fqq 



= 0 ...(Dl) 



When a vector F(m) in the vectors is exchanged by a set of variables (xl, .... 
xq), the plane can be expressed, as explained before, by the following 
expression: 

Ll<m)<xl+ ... + Lq(m)-xq = 0 ...(D2) 
where Ll(m), L(2) Lq(m) are cof actors of each variable (xl. ...,xq). 

By multiplying a normalizing factor N = ((Ll)* + ... -KLq)*)" 1 * to this 
expression and further by multiplying a sign normalizing factor S so that the 
coefficient of one of the variable, for example, xl, to be positive, the following 
expression is obtained. This expression is invariant irrespective to the 
selection of m. 

N-S-Ll(m) xl + ... + NS-Lq(m) *q = 0 _(D3) 



Namely, when the vectors F(l), ... F(q) are in a common plane, the following 
relations stand: 

NSL1(1) = ...=NSLl(q) 

N-S-L2(l) = ... = NSL2(q) 



NSLqd) = ... = N-SLq(q) 

However, when the vectors F(1),...F(2) are not exactly in a common plane, 
but are positioned nearly in a common plane, the following relations stand 

N-S-Ll(l)* ... -NSU(q) 

NSL2(1)^ ... *NSL2(q) 
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NSLq(l)- ... -N-S-Lq(q) 

Hie plane degree index PI is an index indicating in what degree the vectors 

F(l) F(q) can be regard as a plane. The plane degree PI shall be defined so 

that the plane degree index is a function of the direction cosine Lj(p), and the 
value of the index PI approaches to a specified standard value PI(0), (for 
example. PI(0) can be identical to zero), when the proportions of the direction 
cosine distribute closer to each other. Namely the fact that PI is equal to the 
standard PI(0) means that the vectors form a common plane, and the fact 
that PI is nearly equal to the standard PI(0) means that they form nearly a 
common plane 

For example, assuming the mean value of the N-SLk(m) to be <Lk> 
deviation of the N-S *Lk(m) from the mean value to be dLk(m), and the means 
square value of the deviation to be msv, it is preferable to define the plane 
degree index by the means square value msv. 

Further, assuming the maximum of the absolute value of N S Lk(m) - <Lk> 
to be M(m), and the maximum of the M(m) to be M, the plane degree index PI 
can be defined as PI = M. 

PI, PI(0) can be arbitrary defined on the basis of any industrial or 
commercial grounds. 

When the absolute value of the difference I PI-PI(0) I between a plane degree 
index PI obtained from q vectors F(l), .... F(q) and the standard value PI(0) is 
smaller than a predetermined allowance dPI, it can be defined that the 
vectors are nearly in a common plane, and the plane is the sought mapping 
plane. It is possible to output the mean value of the direction cosine 
<Ll>,... f <Lq> as a direction cosine of the vectors F(l),...,F(q). Also the 
allowance dPI can be defined on the basis of any industrial or commercial 
grounds. 

When it is assured that the measuring error of the data d(i) are small, it is 
not always necessary to obtain the plane index PL 
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E. <Cases that cofactors are zero > 

Assuming that the measured data d(i) = (xi, yi, ti) are coordinates of a 
moving body, and (xi, yi) are values of the position coordinate (x, y), and tit 
are values of the time coordinate t Further assuming that the following 
relations stands: 



When a set of base functions are defined as xl = 1, x2 = t, x3 = t 1 , x4 = x, x5 = 
y, for this data d(L), and forming a mapping of these data set d(l) f d(2) v d(3), 
d(4), d(5) into the five dimension vector space [xl, x2. x3, x4, x5], the 
following five vectors can be obtained: 

P(l) = (l.tl. (tl) 1 , xl. yl) 

P(2) = (l,t2, (t2) f . x2, y2) 

P(3)= (1. t3. (t3)«. x3. y3) 

P(4)= (1. t4. (t4) f . x4. y4) 

P(5)= (1. t5. (t5)*. x5. y5) 

Hie cofactors A14. A24, A34. A44, A54, regarding the fourth column and the 
cofactors A15. A25, A35, A45, A55, regarding the fifth column of the 
determinant of the vectors P(l), P(2), P(3), P(4), P(5) axe always equal to zero, 
because of the relations (El) and (E2). 

Hie fact that cofactors regarding a specific column are always equal to zero 
indicates that the selected candidate mathematical relation is not 
appropriate, and some unnecessary base functions are included in the 
mathematical relation. 

As to this example, either of die base functions [yj or [xj is unnecessary For 
example, when a candidate mathematical relation is selected so that the 
data set d(i) can be mapped into a four dimension vector space (l,t,t*,x), a 
relation between the time t and the x co-ordinate can be obtained from the 
direction cosine. And when a candidate mathematical relation is selected so 
that the data set d(i) can be mapped into a four dimension vector space 
[l,t,t*,y1, a relation between the time t and the y co-ordinate can be obtained 
from the direction cosine. 



x + 2t + 3-l=0 
y+4t*+5t + 6-1 = 0 



...(El) 
~.(E2) 
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F. Mathematical basis of the method for obtaining the direction cosine, 
using an eigenvalue> 

Even if the mapping of a set of measured data d(i) through base functions (fl, 
f2, .... £3) could form a plane, when measured data contain errors, and the 
values of base functions also contain errors, thus the mapped paints 
corresponding to the measured values distribute at a region near to a plane. 
In this method, the direction cosine of a plane, the sum of the squares of the 
distance between those points and the plane (i.e., the square of the length of 
the perpendicular lines) is the minimum, is calculated. 

The direction cosine (Ll,...JLq) of a plane in a q dimensional space is 
perpendicular to this plane, and its absolute value is one. Thus the length 

h(i) of the perpendicular line from each point Pi (Fli, F2i Fqi) to this 

plane is given by the inner product <(L1, ... LI) * (Fli, F2i, ... >Fqi)> of the 
direction cosine and the vector Pi. 

h(i) = HFli + ... + Lq-Fqi ... (Fl) 

The sum H of the squares of the lengths h of the perpendicular lines is given 
by the following expression, using a correlation sum <a I b>: 
H =Zh( i ) * 

=£(LlFli + ... + Lq-Fqiy 

q q 

= Z ( Z(ULb-I(Faifl)i))) ...(F2) 
a=l b=l i 
(In this specification, £ is a sum regarding running variable i.) 

Here, a correlation sum is defined as follows: 
Correlation sum < a I b> = HFaiFhi) ...(F3) 

Following expression can be obtained; 

q q 

Zhtt) , = I ( I(La-Lb-<alb>)) ...(F4> 
I a=l b=l 
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Because Ii(l, 2, q) is a direction cosine, the following relation stands: 
(Ll)* + ... + (Lq)* = 1 ...(P5) 



10 



Now, the problem is to find the values of LI, L2, Lq, which give the 
minimum of the square sum of h*, under the condition (F5). This problem can 
be solved using the Lagrange method of undetermined multipliers. In the 
Lagrange method, for example, the maximum or minimum of a two variable 
function f(x, y) under a condition ( g(x, y) - 0 ) can be found, by making a 
estimation function ( F(x, y) = f(x, y) -X-g(x, y) ), and calculating the value of X. 
x, and y, from the conditions (dF/dx = (K dF/ay = 0^ g(x, y)=0). 



15 



The Lagrange method is applied to the present problem, where the variables 

are LI, Lq. The estimation function F(L1 Lq) is defined as follows: 

F (Ll f - . . Lq) =h»- V ((LI) (Lq)*- 1) „.(F6) 



20 



Because the derivative of the estimation function with respect to the 
variables Ll. ...,Lq are equal to zero, the following eigenvalue equation can 
be obtained: 



f<l|l> <1|2> ... <llq* 
<2|1> <2|2> ... <2|q> 

^<q|l> <q!2> ... <qlq>. 




...(F7) 



25 The matrix comprised of correlation sum <Uj> at the left side is called 
correlation sum matrix C, in this Specification. The eigenvalue equation is 
equivalent to the following secular equation: 
<1|1>-X <1|2> ... <l|q> 
<2J1> <2|2>-X... <2|q> 



7 



<q|l> <q|2> 



<qlq>-X 



...(F8) 



When this secular equation can be solved, and assuming that an eigenvalue 
is Xm, and an eigenvector corresponding to the eigenvector is (him, L2m, 
tr, Lqm), the following relation stands: 

H = Lh(i)« 
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w 



20 



25 



30 



40 



^(LlFi+.-. + Lq-Fqi) 1 

q q 

=1 (Z (La-Lb-Z (FaiFbi))) 
5 a=l b=l i 

q q q 

= 1 La- (ILb <a lb>) =ZLaXmLa 
a=l b=l a=l 
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= Xm((Ll) , + ... + (Lq)*) = Xm ... (F9) 



This means that when the square sum of the length of the perpendicular 
lines is minimum, the square sum H of the length of the perpendicular lines 
15 is Xm, when it is zero, the vector points P(i) are on a plane, and when it is 
nearly equal to zero, the vector points are distributed in a region near to a 
plane. When the eigen value Am is nearly equal to zero, the mathematical 
relation between the measured data © can be identified as follows: 

20 LI fl + L2f2 + ... + Lqfq = 0 ...(F10) 



The validity can be estimated whether the eigen value Xm is smaller than a 
predetermined value A©. 

35 

2fi A plurality of the eigenvalues can be simultaneously valid. For example, 
mappings of the measured data distribute on a cross-sectional line of two 
planes, the equations of the two planes corresponding to two eigenvalue can 
be identified as the mathematical relations. 



30 For example, consider a case that the input data d(i) are co-ordinates of a 
45 trajectory (t, x, y), and the data d(i) of the trajectory (t v x, y) are mapped into 

a five dimensional vector space[xl, x2, x3, x4, x5). namely, (t-»xl. t*-»x2, 
l-*x3. x-»x4, y-*x5]. Assuming that the following two mathematical 
relations are obtained from the eigenvalues and the eigenvectors: 
50 .W 2x1+0x2 + 5x3+1x4 + 0x5 = 0 _(F11) 
O x I + 3 x2 + 4«x3 + 0 x4 + 1x5 = 0 ...(F12) 



55 
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Direction cosines can be obtained by multiplying a normalizing factor N t , N 3 , 
respectively to each set of the coefficients (2,0,5, 1,0), and (0,3,4,0,1). 

Ni = (2-2 + 5-5 + M) m 

N,= (3-3 + 4-4+ 1-1 ) 

In this case, the following two mathematical relations can be identified and 
they are simultaneously valid. 

2- t + 5 + x = 0 ... 0F13) 

3- t» + 4+ y =0 ...(F14) 



Sometimes, weighting of measured data is necessary. There is a case that a 
mathematical relation is identified from a time sequential data d(i), 
imposing a light weighting to old data and a heavy weighting to new data. 
For example, when a temperature shall be controlled to follow a 
predetermined curve by controlling the heat supply, there is a case that it is 
not preferable to treat old and recent data with an identical weighting. 

As explained before, the length hQ) of the perpendicular line from each point 
Pi(Fli F2i, ....Fqi) to a plane can be expressed by an inner product of the 
direction cosine and the vector Pi: 

h© - Ll-Fli + ... + Lq Fqi ~. (Fl) 



Multiplying a weight Wi to the value h t following expression can be obtained: 
Wi-hCi) = Wi(LlFli + ... + Lq-Fcp) ... (F15) 

The square sum H of the lengths h(i) of the perpendicular lines with the 
weighting is given by the following expression: 
H = £Wi 2 h(i)* 
= IWi* (Ll Fli + ... + Lq-Fqi) 2 
q q 

= Z < £ (U Lb-EWi*- (Fai-Fbi))) ... (FIG) 

a=l b=l i 

Here, the definition of the correlation sum is expanded as follows: 

Correlation sum <a I b> = £Wi 2 . (Fai-Fbi) ... (Fl 7) 

( The definition by (F3) corresponds to a case that all the Wi are equal to 
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one.) 

The following relation stands: 
H=IWiMiG)* 

q q 

=1 ( I ( LaLb-<alb>)) ...(F18) 
a=l b=l 

The correlation sum defined by (F17) is a real symmetrical matrix. And the 
expression (F18) has same form as the expression (F4), therefore, all of the 
arguments explained on the ground of the definition of the correlation sum 
according to (F2) are valid for this definition. 

G. <The eigenvalues of the correlation sum matrix are positive real 
numhers> 

1) The correlation sum matrix C is a real symmetric matrix. Therefore the 
matrix is an Hermitian matrix, and all of its eigenvalue are real number. 

2) It is known that the necessary and sufficient condition for that all of the 
eigenvalues of an Hermitian matrix T are positive is that for any non-zero 
vector x, the scalar product <xT-x T > is non-negative. 

According to the definition of a correlation sum matrix (F3): 

<a I b> =Z (Fai-Fbi) ... (F3> 

The scalar product is obtained, when a vector (xl, x2 xq) and its 

transposed matrix (xl, x2, xq) T are multiplied to a correlation sum matrix 
C, respectively from the left side and the right side. Using the expression 
(F9), the scalar product <xC-x T > can be transformed as follows: 
<xCx r > = ££xa-xb^a|b> 

q q 

= Z (ZOta xb-HFai Fbi))) 
a=l b=l i 

^ I (xl-Fli + ... + xqFqi)'£0 
i 
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As shown here, <xC-x T > is non-negative. Also when the correlation sum 
matrix is defined by (F17). it can be proved that the all the eigenvalues are 
non-negative real numbers, in the same manner. 

H. <A method for obtaining the eigenvalues between a predetermined 
value Xo and the zero> 

Because all the eigenvalues of a correlation sum matrix are non-negative, all 
the eigenvalue smaller than the predetermined value X<, are between the zero 
and the value A* An eigenvalue in this region can be found as follows, for 
example. 

Assuming that a matrix K having q rows and q columns is given, as follows; 
'Kll K12 ... Klq\ 




K21 K22 ... K2q 
Kq2 ... Kqq i 



(HI) 



Selecting n diagonal elements Kii(i=nl, n2, nn) from this matrix, a 
determinant K(nl, n2, nn) can be formed, which is composed of the nl-th, 
n2-th v nn-th columns and rows of the matrix K. As an example, matrix 
K(l,3) is shown below: 

K(1,3)=IK11 K13I ~-(H2) 
JK31 K33| 

A trace(n) is defined as a sum of all the determinant K(n 1, n2, .... nn), which 
include n diagonal elements of the matrix K. Here, trace(O) is defined to be 1. 

trace(0) = 1 ... (H3a) 

trace(l) = Kll + K22 + ... + Kqq ... (H3b) 

trace(2) = I IKiiKijl 

i*j|KjiKjj ...<H3c) 
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Kii 


Kii 


Kik 




trace<3) = 


i I 


Kji 


Kjj 


Kile 






Kki 




tin i 

Kkk 














k*i 












Kll 




K12 




Kiq 




K21 




K22 




K2q 


trace(q) = 


Kql 




Kq2 




Kqq 



(H3d) 



(H3e) 



The secular equation corresponding to the matrix K can be expanded as a 
polynomial of the powers of (-X), as follows: 



0 = 



Kll-X. K12 
K21 K22 X 



Kq2 



Klq 
K2q 

Kqq-X 



...(H4) 



Kql 
Q 

= Z(trace(q-n))((A)-) 
n = 0 

This is a polynomial of (-X). Therefore, it is easy to find a zero point in a 
region between zero and \ using a conventional method, for example, by 
dividing the region into small sections. 



I. < An application of the identification of the mathematical relation to an 
analysis of a trajectory of a golf head> 

Quantitative diagnosis of golf swing is possible. For example, finding die 
direction cosines of planes, in which the golf club head moves in the initial 
period, middle period and the final period of the golf swing, one can diagnose 
whether the swing trajectory is in a plane or not, according to the change of 
the direction cosines in the swing, and diagnose how the club head is 
accelerated. 



For example, assuming that a trajectory of the golf club head is an 
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intersectional curve of a fixed plane and an ellipsoid, which has a fixed 
center and its size changes with time t according to a second order function 
of time t, the trajectory can be expressed as follows: 

#l)p-x* + qy , + r-z t + a-x + by + c-z 
+d0 + dlt + d2t* =0 ...(ID 

#2) a x + b y + c-z + d-l ■ 0 ...(12) 

It is possible to determine the values of the parameters p, q, r, a, b, c, dO, dl, 
d2 9 dO of the expressions (II) and (12), by selecting an appropriate base 
functions, a candidate mathematical relation and by mapping into a vector 
space F, as explained in connection with the expressions (6) and (7). 

A method to find the change of the angular position, namely, the angular 
velocity and the angular acceleration of this movement in this plane (a x + 
by + c-z + d « 0) is explained below: 

Here, a vector P(xi, yi, zi) is a vector in a real vector space (x. y. z ), which is 
not an imaginary vector space F. Assuming that the a second plane is a plane 
parallel to this plane (a x + by + c-z + d = 0) and the origin of this real vector 
space is included in the second plane. The distance D(i) between the 
measured position of the club head and the second plane can be calculated 
from the measured valued d(i). D(i) is a scalar quantity. The direction cosine 
N of the second plane in the real vector space is proportional to (a, b, c), 
namely N = k(a, b, c). Where k is a normalizing factor and is expressed as 
follows: 

k = (a t + b > + c t ) •« ...03) 

A vector Q of a point, which is an intersection of the second plane (a x + b y + 
c-z + d =0) and a normal line of the second plane passing through the vector 
" can be expressed as follows: 

Q(i) = P(i)-DC0-N .-.(14) 

Setting a reference vector Q, for example Q(l), and calculating an inner 
product (ie., scalar product) <Q(1) -Q(i)> of the vectors Q(l) and Q(i). the 
angle e(i) between those vectors can be calculated from the proportion r of 
the product of the absolute value of these vectors IQ(l)l !Q(')I and the 
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inner product: 

r«<QO)-Q6WIQO)hlQC9l 

6(i)=cos- l r ...05) 

When the following mathematical model is supposed, 

8 = a + bt + ct* ...06) 
the coefficients a, b, c can be determined, selecting a candidate mathematical 
relation, which is composed of the following base functions gl, g2, g3, and g4, 
and applying the method according to the present invention: 
gl = 6 G). g2=l, g3=t g4 = t l 

When the coefficients a, b, c of (16) are determined, the angular velocity and 
the angular acceleration can be calculated. 

Next, we consider a third model #3). in which the center of the ellipsoid 
moves with time t according to a first order function of time t, and also the 
size changes with time t according to a first order function of time t, and a 
fourth model #4), in which the direction cosine (a, b, c) changes with time t 
according to a first order function, and also the distance g changes with time 
t according to a first order function of time t. 
#3)p(x-a0 + aO + alt) 1 + q(y .b0 + b2t) f 

+ r•(^-cO + c^t) , (rO+rl-t)* = 0 ... (17) 

#4)(a + bt)-x + (c + dt)y + (e + ft)-z-(g + ht) =0 ... (18) 

The parameters of the model #3) can be determined by selecting a candidate 
mathematical relation, corresponding to the following base functions: 
fl»*,fi = y t l ffi = s' > f4 = t\ 

f5 = x. f6=y r f7 = z, fS = t, 

f9 = xt, flO = yt, fll=zt. fl2 = l 

The parameters of the model #4) can be determined by selecting a candidate 
mathematical relation, corresponding to the following base functions: 
fl = x, f2 = y, C3 = z, f4 = t, 

f5 = t x, f6-tt, t7 = t-z, f8=l 

J. < An example of identification of a differential equation from measured 
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data> 

Assuming that it is known that the measure data d = (x, y) can expressed by 
the following expression: 

a-dy/dx + by + c-y* +dx + e-x* +f-l= 0 ... (Jl) 

The parameters a, b, c f d, e, f can be determined by selecting following base 
functions fl, £2, .... f6: 

fl = dy/dx = (y(i) - y(i-l))/(t(i)-t(i-l)) 

f2 = y = y(i) 

f3 = y*=y(i)* 
f4 = x = x(i) 
f5 = x* =*<i)' 
16=1 

The candidate mathematical relation is: 
a-fl + b£2 + c-f3 + df4 + e-f5 + f-ffi = 0 

The parameters a, b, c. d, e t f can be determined from the direction cosine in 
a vector space spanned by these base functions. However this method is 
applicable only when the density of data is sufficiently thick for neglecting 
the errors caused by the approximation of the derivative by a division. 



Claims 

1 . An apparatus for identifying a mathematical relation between base variables (x1 1 x2 xp), when a set of input data 

d is composed of p base variables (x1 , x2 xp), and a plurality of data set d(i) of such a data set d is inputted, the 

data sets d(i) are distinguished by an input data distinguishing parameter V, comprising: 

a data acquisition means (1) for acquiring values xji of base variables xj (j = 1,...,p) in sets of input data d(i) (i 
= 1,2 n), which are distinguished by the input data distinguishing parameter i; 

a base function defining means (4) for selecting a mathematical function program (m = k ), which is one of the 
mathematical function programs gm distinguished by a mathematical function specifying parameter m, the 
mathematical function programs output a function value corresponding to their input references, the base func- 
tion defining means (4) defines a base functions fk, by specifying a base variable to each of one or more than 
one input references of the selected mathematical function program fk; 

a candidate mathematical relation specifying means (5), which specifies a candidate mathematical relation, by 

specifying a plurality of base functions fk (k = 1, 2 q) from the base functions defined by the base function 

defining means; 

a vector component acquisition means (6) for sending the values of the base variables acquired by data acqui- 
sition means (1), to each of the input references of base functions fk (k = 1, ... q) included in the candidate 
mathematical relation specifying means, which is defined by the candidate mathematical relation, for each data 
set d(i), and for acquiring and storing the output value Fki of the base functions fk into an array VC_ARRAY t in 
which the values of the functions Fki are stored in a form of a vector F(i) in a q dimensional space, for each 
input data specifying parameter V. 

a direction cosine acquisition means (8) for acquiring the direction cosine (L1 Lq) of a mapping plane 

spanned by the points P(i) in a vector space which correspond to the vectors stored in the vector component 
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array VC_ARRAY t wherein the direction cosine acquisition means (8) comprises: 

a correlation sum calculating means (8B1) for calculating the correlation sum (alb) of the vector compo- 
nent array VE_ARRAY; 

5 an eignevalue calculating means (8B2) for calculating the minimum eigenvalue among the eigenvalues X 
(1),..., X (q) of a correlation sum matrix, the elements of which are correlation sums (alb) ; 
and an eigenvector calculating means (8B3), which calculates the eigenvector corresponding to the mini- 
mum eignevalue, and outputs the eigenvector as the direction cosine (L1 Lq); 

w and a mathematical relation outputting means (10) for outputting a mathematical relation corresponding to the 

following expression or an mathematical expression or values, which can be deduced from the following 
expression, when the direction cosine of the plane corresponding to each base function is Lk( k = 1 , 2, q): 

q 

75 £(Lk-fk) = 0 

k=1 



2. An apparatus for identifying a mathematical relation according to claim 1 , wherein the eigenvalue calculating 
20 means (8B2) multiplies a weighting Wi 2 to the product of elements Fai and Fbi in the may VC_ARRAY f then obtains 
the sum of them in respect with the parameter i to obtain the correlation sum < alb ) . 

felb)=ZWi 2 *(Fai'Fbi) 

25 3. An apparatus for identifying a mathematical relation according to claim 1 , wherein further comprising: 

a validity estimating means (9) for judging whether the obtained eigenvalue is larger than a predetermined 
value or not; 

and a file (F1 1) storing a control program which controls so that when the validity estimating means (9) does 
not estimate that the eigenvalue is smaller than the predetermined value, the candidate mathematical relation 
specifying means (5) selects another candidate mathematical relation, and when the validity estimating means 
(9) estimates that the eigenvalue is smaller than the predetermined value, and the direction cosine correspond- 
ing to each of the base functions is Lk, the mathematical relation outputting means outputs a mathematical 
relation corresponding to the following mathematical relation of deduced from the mathematical relation. 

q 

£(Lk-fk) = 0 

k=1 



4. An apparatus for identifying a mathematical relation according to claim 1 , wherein the apparatus is incorporated in 
a measuring apparatus, and used as an apparatus for interpolating or extrapolating data. 

5. An apparatus for identifying a mathematical relation according to claim 1 , wherein the apparatus is incorporated in 
45 an apparatus for identifying trajectories of a moving object, and acquires input data of position as a function of a 

time parameter Y, then outputs one or a plurality of mathematical relation(s), for identifying the trajectory of the 
moving object. 

6. An apparatus for identifying a mathematical relation according to claim 1 , wherein the apparatus is incorporated in 
so an apparatus for identifying a form of a curve or a surface, and acquires input data of coordinates of infinite points 

on the curve or the surface, then outputs one or a plurality of mathematical relation(s), for identifying the form of the 
curve or the surface. 

7. An apparatus for identifying a mathematical relation according to claim 6 f wherein the form of curve or surface is a 
55 form of an object, and coordinates of points on the object are given as input data d(i). 

8. An apparatus for identifying a mathematical relation according to claim 6, wherein the form of curve or surface is a 
form of curve or surface in an abstract space, which represents a mathematical relation of variables, and values of 
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the variables are given as input data d(i). 

A memory media, which can be read out by a computer, and is used in an apparatus for identifying a mathematical 
relation according to claim 1 , the memory media contains the following files: 

a file (F1) storing a data acquisition program for acquiring sets of input data d(i); 

- a file (F3) storing a plurality of mathematical function program, which can output a function value correspond- 
ing to a set of its input variables: 

a file (F4) storing a base function defining program for defining base functions and/or defined base functions; 

- a file (F5) storing a candidate mathematical relation specifying program for specifying a candidate mathemat- 
ical relation and/or a set of base functions as a specified candidate mathematical relation; 

a file (F6) storing a vector component acquisition program, which sends the values of base variables to the 
input references of the base functions fk (k = 1 ,...,q) included in the candidate mathematical relation specified 
by the candidate mathematical relation specifying program, for each input data specifying parameter i, and 
acquires the value Fki of the base function to store a vector F(i), the component of which is the function values 
Fki (k = 1, ...,q), into a vector component array VC_ARRAY, for each input data specifying parameter i; 
a file (F8) storing a direction cosine acquiring program, which acquires the direction cosine of a mapping plane 
corresponding to the vectors F(i) in the vector component array VC__ARRAY; 

- and a file (F10) storing a mathematical relation outputting program for outputting a mathematical relation, 
which corresponds to the following expression, when the direction cosine corresponding to each of the base 
functions fk is Lk (k = 1 , .... q): 

£(Lk.fk) = 0 

k=1 
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Fig. 3 
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Fig. 5 
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Fig .7 

d (1) = ( Pi.Ti.Xi Yi Zi) 
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